Measurement of intrinsic Dirac fermion cooling on the surface of a topological 
insulator Bi 2 Se 3 using time- and angle-resolved photoemission spectroscopy 
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We perform time- and angle-resolved photoemission spectroscopy of a prototypical topological 
insulator Bi2Se3 to study the ultrafast dynamics of surface and bulk electrons after photo-excitation. 
By analyzing the evolution of surface states and bulk band spectra, we obtain their electronic 
temperature and chemical potential relaxation dynamics separately. These dynamics reveal strong 
phonon- assisted surface-bulk coupling at high lattice temperature and total suppression of inelastic 
scattering between the surface and the bulk at low lattice temperature. In this low temperature 
regime, the unique cooling of Dirac fermions in TI by acoustic phonons is manifested through a 
power law dependence of the surface temperature decay rate on carrier density. 



Electrons on the surface of a three-dimensional topo- 
logical insulator (TI) are massless Dirac fermions with 
linear energy- momentum dispersion These elec- 

tronic systems are promising for novel applications rang- 
ing from spin-based field effect transistors Q, ultrafast 
opto-spintronic devices [5[ and high-speed quantum infor- 
mation processors [6], whose performance depends crit- 
ically on the dynamics of hot carriers. Unlike the case 
in graphene [7[, helical Dirac fermions in a TI interact 
not only with phonons but also with an underlying bulk 
reservoir of electrons [3[ . Therefore it is important to un- 
derstand separately their coupling mechanisms to each of 
these degrees of freedom. However, high frequency op- 
tical conductivity [8| and ultrafast optical experiments 
do not directly separate bulk and surface signals 
nor do they distinguish different relaxation channels. 

Time- and angle-resolved photoemission spectroscopy 
(TrARPES) is a powerful technique to study carrier dy- 
namics with energy and momentum resolution. A recent 
TrARPES study has shown that persistent surface carrier 
population in a TI can be induced by photo-excitation 
[111 ] . In this letter, we focus on disentangling different re- 
laxation channels following photo-excitation over a range 
of electron doping suitable for device applications. Anal- 
ysis of the TrARPES spectra enables us to separately 
obtain temperature and chemical potential relaxation of 
both surface and bulk. Our data reveal that interband in- 
elastic electron-electron (e-e) scattering is suppressed and 
surface-to-bulk coupling is mediated entirely by phonons. 
At low lattice temperature, this latter scattering channel 
turns off and we observe that the surface temperature 
decay rate scales as a power law with doping indicating 
acoustic phonon assisted cooling of 2D Dirac fermions. 

Single crystals of Bi2Se3 are cleaved under ultra-high- 
vacuum (< 6 x 10 -11 Torr) [12]. Their doping levels 
are determined by measuring their static ARPES spec- 
tra [l2| . TrARPES is performed following a pump-probe 
scheme [Hj that uses 1.55 eV pulses with a fluence of 
100 /iJ/cm 2 as a pump and 6.2 eV ultrafast laser pulses 
as a probe. A time-of- flight electron energy analyzer [14J 



is used to simultaneously collect the energy-momentum 
spectra I(E,k x ,k y ,t) of Bi2Se3 at variable pump-probe 
time delay t without sample or detector rotation [15j. 

To understand the relaxation dynamics of the states 
near the Fermi level Ep, we show low energy ARPES 
spectra along the V — M direction of the surface Bril- 
louin zone at several time delays [Figs. l(a)-(d)]. Fig. 
1(a) shows the surface states (SS) and conduction bulk 
(CB) band of a Bi2Se3 sample with the Dirac point Ejj = 

0. 28 eV below Ep before the pump pulse arrives (t = - 
5 ps). The measured bandstructure is consistent with 
prior studies [16]. The ARPES spectra taken immedi- 
ately after the pump excitation [Fig. 1(b)] resemble a 
thermalized electron distribution at elevated temperature 
for both SS and CB. This hot distribution cools down 
progressively at longer time delays [Figs. 1(c) and (d)]. 
However, it can be seen from the difference of momentum 
integrated spectra AI(E,i) = I{E,t) - I(E,t = —bps) 
that the equilibrium has not been reached within the 
measurement time window [Fig. 1(e)]. 

In order to quantitatively describe the relaxation dy- 
namics of the hot bulk and surface populations in Fig. 

1, we separately analyze the spectral intensities for the 
surface and the bulk. Fig. [2^a) shows I(E, t) for SS and 
CB at t = 0.25 ps obtained in their respective energy- 
momentum region [Fig. 1(a)]. We fit I{E,t) for both 
SS and CB with the following equation [17j to extract 
their respective electronic temperatures T e (t) and chem- 
ical potentials fi(t) as a function of time: 

/oo 
de[f FD (E,T e (t),n(t))xD(E)] 
-CO 

G(E-e,w(t)) (1) 

Here A(t) is a scaling factor |17j , / FD is the Fermi-Dirac 
distribution, G(E, w) is a Gaussian function with energy 
independent width w that includes the instrumental res- 
olution and time-dependent broadening due to increased 
scattering rate after photo-excitation [18|. The time in- 
dependent D, which is a product of density of states 
(DOS) and photoemission matrix element [19|, is ob- 
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FIG. 1. TrARPES spectra of Bi2Se3. (a)-(d) Energy- momentum spectra sliced through 3D data volume along T-M direction at 
-5 ps (a), 0.25 ps (b), 2.5 ps (c) and 15 ps (d) taken at 15 K on an Ed = 0.28 eV sample, (e) Momentum integrated difference 
spectra subtracting the spectrum at t = -5 ps. The blue dashed lines in (a) define the boundary between conduction bulk band 
(CB) and surface state (SS) for the analysis performed in the text. 
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FIG. 2. Fitting of momentum-integrated spectra (a) Momen- 
tum integrated spectra for SS (green) and CB (red) and their 
respective fits (purple and orange lines) at t = 0.25 ps for 
the data shown in Fig. 1. (b) Constant energy slices of the 
difference spectra as a function of t at 0.1 eV for the same 
data and fit. The spectra are normalized over the pre-pump 
intensity Iq = I(t — — 5 ps). 



tained through a global best fit for all t [12]. Both the SS 
and CB fittings agree very well with data for all t [Fig. 
[2](b)] (also see [HI). The fact that both SS and CB pop- 
ulations can be described by Fermi-Dirac distributions at 
t =0.25 ps [Fig. Efa)] suggests that intraband thermal- 
ization of SS and CB electrons is established within the 
instrumental resolution time (~ 200 fs). 

Fig. [3] shows the obtained electronic temperatures and 
chemical potentials for both the SS and CB. We focus on 
the temperature first. As shown in Figs. [3ja), SS (CB) 
temperature T e ss (T e CB ) starts at the lattice temperature 
Tl = 15 K at t = -5 ps. The maximum temperature 
increase of ~ HOOK at t = ps allows us to estimate 
an electron specific- heat coefficient 7 of 2 mJ/rao/ • K 2 
(assuming 50 nm penetration depth of the pump pulse 
[Hi). The relaxation time agrees with previous optical 
studies [To| . 



The analysis of the time dependence of the chemi- 
cal potential reveals carrier population dynamics both 
within and between the bands. In general, chemical po- 
tential is a function of both temperature and particle 
number N [2cj. In the case where N is fixed, we de- 
rive [12] a relationship between chemical potential and 
our measured temperature assuming linear and parabolic 
dispersion for the SS and CB bands [20] respectively. 
The SS and CB chemical potentials derived this way 
(fj^ and /ijy ) are plotted in Fig. [3jb) along with the 
chemical potentials directly obtained by fitting our data 
to Eq. (1) (/i ss and /i CB ). Remarkably, a good overlap 
between /i ss and /i^ as a function of t is observed at 
Tl = 15 K [Figs. [3fb) green] indicating constant SS 
particle number N ss . ji CB is significantly higher than 
/ijy [Figs. [3f b) red], suggesting instantaneous 7V CB in- 
crease upon photo-excitation. The increasing N CB shows 
that photo-excitation from deeper valence bands primar- 
ily populates CB, consistent with a related st udy that 
shows bulk dominant direct optical transition [11]. The 
fact that N ss stays constant indicates that there is no net 
particle transfer from CB to SS. 

In order to understand whether the transfer of carriers 
between surface and bulk can happen at higher lattice 
temperature, we now investigate the chemical potentials 
and temperature dynamics at 300 K [Figs, [3fc) and (d)]. 
The most striking feature is that ja ss is much bigger than 
jn s ^ [Fig. [3jd)], indicating a strong increase of surface 
carrier density in contrast to 15 K [Fig. [3fb)] . To find 
where the extra surface carriers come from, we notice 
/i CB is lower at 300 K than at 15 K [Fig. [3{b) and (d) 
solid red traces]. By calculating the change of particle 
numbers after photo-excitation of SS and CB (AN SS and 
AN CB ) [H[, we see that the increase of AiV ss at 300 K has 
a slower rise time than AN CB at 15 K [Fig. [3](d) inset]. 
This shows that SS particle number increase is not due 
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FIG. 3. Lattice temperature dependent interband surface-bulk electron dynamics, (a) Electronic temperature T e and (b) 
chemical potential fi for SS (green) and CB (red) at lattice temperature Tl = 15 K. (c) and (d) are at Tl = 300 K. The green 
dashed traces are the SS chemical potential as a function of T| s under constant particle number. The red dashed traces 
are /i^ for CB. The change of particle number is shown in (d) inset. (The features before t — are due to the processes where 
the roles of pump and probe are reversed, which is commonly seen in TrARPES [131].) 



to direct optical transition. Rather it is a result of trans- 
fer from CB to SS through phonon scattering at elevated 
lattice temperature due to higher phonon-scattering rate 
above Debye temperature (Q D = 182 K (HI) [H|. The 
evidence for such phonon-assisted CB-SS scattering can 
also be seen in the temperature dynamics [Figs. [3ja) 
and (c)]. Tf and T e CB only equilibrate at 300 K within 
the probed time window, indicating that phonons scatter 
carriers, transferring energy and mediating the thermal- 
ization between SS and CB at 300 K. Since inelastic e-e 
scattering is mainly responsible for ultrafast thermaliza- 
tion [22| within a single band, the lack of thermalization 
between SS and CB at 15 K [Figs. [3ja)] suggests that 
inelastic e-e scattering time between SS and CB exceeds 
the measured time window. (See [12] for more T e ss vs. 
T e CB at 15 K.) Such suppression of interband inelastic e-e 
scattering is likely due to the kinematics constraints [23| . 

Having shown that interactions between SS and CB are 
suppressed at 15 K [Figs. [3fa) and (b)], we proceed to in- 
vestigate their respective intraband cooling mechanisms. 



Because different mechanisms can be distinguished based 
on carrier density dependence, as demonstrated in other 
Dirac materials [24|, we measure the intraband cooling 



dynamics of Bi2Se3 over a wide range of dopings [12 ] . 
To vary the doping, we used samples from two different 
batches with high and low carrier concentrations along 



with surface doping effect 25, 26]. The surface doping 
effect has been commonly observed in topological insu- 
lators such as Bi2Se3 [25] and is attributed to either ad- 
sorption or moving of Se vacancies towards the surface 
[26|. Figure ffl[a) shows that the surface temperature has 
two decay components r^ s and r| s . The fast-decaying 
component r^ s is independent of Ejj. This is consistent 
with optical phonon cooling, which is insensitive to dop- 
ing as observed in graphene [l?} . When the electronic 
temperature cools below ~ 600K, the optical phonon 
cooling becomes less effective. Instead, the decay be- 
comes dominated by a density dependent component r| s 
whose rate increases with electron doping [Fig. Ufa)]. In 
contrast, Figure H^b) shows that bulk temperature dy- 
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FIG. 4. Electron density dependent intraband cooling of sur- 
face Dirac fermions. (a) Electronic temperature as a function 
of t of surface electrons with different Ed at 15 K. Different 
colors corresponds to Ed in matching color in the inset, (b) 
T CB at 15 K. Inset: Inverse cooling time of the slow compo- 
nent of T e as a function of Ed at 15 K for SS (green triangles) 
and CB (red squares). l/r| s oc Ed 3 at 15 K. 



namics is qualitatively different. The fast component 
rf B = 0.7±0.1 ps [Fig. H^b)] matches with a previous 
ultrafast study, which has been attributed to an opti- 
cal phonon mediated intraband cooling of CB [10]. The 
slow component r§ B exhibits no discernable dependence 
onE D [Fig. lib)]. 

We use the density dependence of the slow component 
to reveal its microscopic origin. We extract r| s through 
a biexponential fit of Tf(t) [Hi and plot it as a function 
of doping Ed in Fig. HJb) inset. Such exponential decay 
of the electronic temperature in a Dirac system has been 
previously predicted for graphene [2J a nd also observed 
in graphite [HI and recently Bi2Se3 [ll[. The SS cooling 
rate scales as l/r^ oc i?£> 3±0 ' 5 , which is in good agree- 
ment with the theoretical calculations of cooling of 2D 
Dirac fermions via acoustic phonons [24j. In stark con- 
trast, 1/t2 B exhibits little dependence on doping, which 
agrees with theoretical predictions that acoustic phonon- 
mediated carrier cooling in 3D metals is density indepen- 
dent [29]. Consistent with the phonon assisted surface 
bulk coupling presented earlier, the cooling rate at 300 



K lies in between that of the SS and CB at 15 K [l2j|. It 
is worth mentioning that in the low doping regime where 
fi « ksT e , acoustic phonon cooling of Dirac fermions 
will be further suppressed and is predicted to have a dif- 
ferent decay dynamics {24j|. As a result, the doping de- 
pendence of the slow decay component may deviate from 
the power at low doping levels. 

In conclusion, we have directly visualized different 
scattering channels of surface and bulk electrons on a 
TI with TrARPES. At high temperature, phonons are 
largely responsible for scattering electrons between the 
surface states and the conduction bulk band. At low 
temperature, the coupling between surface and bulk is 
suppressed. We further reveal the surface cooling rate 
follows a power law dependence on the carrier density, 
which is a signature of acoustic phonon mediated cool- 
ing of 2D Dirac fermions. This suggests that by tuning 
to the charge neutrality point one can enter a regime 
where acoustic phonon scattering of Dirac fermions is 
completely eliminated and transport properties are deter- 
mined solely by structural and chemical defects 
An exciting possibility is the creation of long-lived hot 
photo-carriers that can be used for high-efficiency pho- 
tothermoelectric applications [H|, [HI . 
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